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INTRCDUCTION 

Computer  codes  for  calculating  realistic  detonation  properties  are 
important  for  predicting  the  performance  of  high  explosives  and  for  as- 
sessing the  usefulness  of  new  explosive  formulations.  Furthermore, 
tremendous  cost  and  manpower  savings  can  be  realized  by  the  judicious  use 
of  computer  codes  employing  an  equation  of  state  that  realistically  pre- 
dicts detonation  properties.  By  properly  employing  the  carputer  model, 
the  chemist  can  eliminate  or  reduce  research  areas  shewing  little  promise 
(thus  saving  synthesis  and  property  evaluation  costs)  and  concentrate  on 
research  areas  offering  the  greatest  chance  of  success. 

The  most  widely  used  equation  of  state  for  predicting  the  detonation 
properties  of  condensed  explosives  is  the  semienpirical  Becker-Kistiakcwsky- 
Wilson  (BKW)  equation  of  state.  This  equation  of  state  has  been  incor- 
porated in  thermo-hydro  codes  such  as  Fortran-BKw\  RLJBY^,  and  TIGER^. 

4 

The  original  equation,  based  on  a Kistiakowsky  and  Wilson  modification 

567  89  10 

of  an  equation  developed  by  Becker  ' ' , has  undergone  many  adaptions  ' ' , 

before  obtaining  its  current,  widely  used  formP 

S'1  + X*SX  (1) 


where  x = < £ nJ^/[V(TM-9)  ] 

In  equation  (1)  P is  the  pressure  of  the  system,  V is  the  volume,  T is 
the  absolute  tenperature,  and  R is  the  ideal  gas  constant.  The  parameters 
a,  8,  x,  and  0 are  etpirical  in  nature  while  k^  is  the  covolume  factor 
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of  the  i detonation  product  species  and  m is  the  nunber  of  moles  of 

the  i^  detonation  product  species. 

This  equation  was  initially  fit  to  experimental  detonation  data  by 
11  12 

Ccvan  and  Fickett  in  1956.  later  Mader  re-parameterized  BKW  to  give 
the  set  of  parameters  (a,  3,  k,  0,  k^'s)  presently  used  in  the  Fortran- 
BKW  code.  These  are  shown  in  Thbles  I and  II.  Slightly  different 
parameters  are  incorporated  in  the  RUBY  and  TIGER  codes. 

To  determine  the  best  BKW  parameters,  one  adjusts  a,  B,  < and  6 and 
the  covolume  factors,  k^,  of  the  appropriate  detonation  product  species 
to  fit  the  best  explosive  performance  data  available.  An  estimate  of 

the  covolume  factor  of  each  species  is  obtainable  through  the  hard  sphere 

* 

repulsion  potential,  r , for  that  species.  Hie  covolume  factor  and  hard 

sphere  repulsion  potential  are  geometrically  related.  In  the  early  1960s, 
12 

Mader  adjusted  the  BKW  parameters  to  fit  five  experimental  measurements 
involving  two  explosives,  trinitrotoluene  (1WT)  and  cyclotrimethylene 
trinitramine  (RDX) , and  experimental  Hugoniot  data  for  the  detonation 
species  H^O,  N^,  and  CO^ . Using  this  limited  data  set,  Mader  found  that 
two  sets  of  parameters  (o,  $,  k,  0)  were  required.  (See  Table  I.)  One 
set  of  parameters  was  needed  for  high  density  explosives  whose  detonation 
products  contain  a large  amount  of  solid  carbon  (TNT- type) . A second 
set  of  parameters  was  required  for  explosives  whose  detonation  products 
contain  a relatively  small  amount  of  solid  carbon  (RDX-type) . Using 
these  two  parameter  sets,  the  predictions  of  relative  explosive  perform- 
ance for  C,H,N,0  explosives  have  been  fairly  reliable.  Furthermore, 
absolute  detonation  pressures  and  velocities  are  calculated  reasonably 
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well  as  long  as  the  explosive  is  a C,H,N,0-type,  is  ideal,  is  near  or 
slightly  lewer  than  00  balance  and  is  at  a high  percentage  of  the 
theoretical  maximum  density.  However,  on  an  absolute  scale  the  BKW 
predictions  have  not  generally  been  accurate.  Furthermore,  where 
extremes  in  explosive  c expositions  are  investigated,  predictions  of  the 
relative  behavior  are  also  in  error.  This  error  in  relative  behavior 

1 

is  also  found  when  BKW  predictions  are  made  for  non-C,H,N,0  explosives. 

The  greatest  deficiency  of  the  BKW  Equation  of  State  is  that  it  con- 
tains no  attraction  potential.  For  this  reason  it  can  be  calibrated 

to  a single  point  in  PVT-space  (as  a CJ  point) ; however,  it  then  fits 

g: 

*? 

badly  at  all  other  points.  Furthermore,-  the  farther  the  calculated  point 
is  frem  the  point  of  calibration,  the  worse  the  fit  will  be.  Thus  it  is 
impossible  to  predict  isentrope  shapes  with  the  BKW  Equation  of  State. 
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Table  I.  Fortran-BWK  Parameters  for  Explosives3 


BKW  Parameter 

TNT-type  Value3 

RDX-type  Value3 

a 

0.50 

0.50 

B 

0.09585 

0.16  - 

K 

12.685 

10.91 

e 

400. 

400. 

3 Reference  12. 


A second  equation  of  state  for  predicting  detonation  properties  is 

the  Jacobs -Ccwperthwaite-Zwisler  (JCZ)  equation  of  state.  The  equation 

of  state,  developed  from  first  principles  by  Dr.  Sigmund  Jacobs  , has 

recently  been  incorporated  into  the  TIGER  computer  code  by  Cowperthwaite 
14  15 

and  Zwisler  ' . It  is  anticipated  that  the  JCZ  equation  of  state,  when 

Drooerlv  parameterized,  will  offer  vast  improvements  over  the  semi-empiri- 
cal BKW  equation  of  state  in  the  calculation  of  not  only  the  pressure 
and  velocity  of  detonation,  but  also  the  detonation  temperature,  product 
composition  and  isentropic  expansion  from  the  detonation  state.  Unlike 
the  BKW  equation  of  state,  the  JCZ  equation  of  state  incorporates  both  a 
tempera ture- independent  pressure  term  and  a pressure  term  resulting  from 
thermal  motion.  The  generalized  JCZ  equation  of  state  can  be  represented 
as 

P = P (V,n. n ) + G(T,V,n, n )NRT/V  (2) 

O 1 s s 

where  P is  the  total  pressure  of  the  system,  Pq  is  the  tenperature- 
independent  pressure  term,  V is  the  volume,  m’s  are  the  moles  of  the 
individual  product  species,  G is  the  Griineisen  parameter,  R is  the  ideal 
gas  constant,  N is  the  total  moles  of  product  species,  and  T is  the 
absolute  temperature  of  the  system.  Furthermore,  the  Griineisen  parameter 
G is  composed  of  several  terms  - one  of  which  is  a potential  function. 

The  JCZ2  equation  of  state  incorporates  the  MIE  potential  function  while 
the  JCZ3  equation  of  state,  the  subject  of  this  report,  employs  an 
expo-p,m  potential  function.  In  the  expo-p,m  potential  function 


f, 

i 

f I 

■ 

or 

if  = <t>r  + <J>a  (4) 

<t>  is  the  total  potential  for  the  pair  of  molecules,  <Pr  is  the  repulsion 
term,  and  <p  is  the  attraction  term.  In  equation  (3)  £ is  the  absolute 

u 

value  of  the  potential  at  its  minimum,  r is  the  molecular  separation  at 

★ 

the  specified  potential,  r is  the  molecular  separation  at  the  potential 

minimum,  p is  the  repulsion  parameter,  and  m is  the  attraction  parameter. 

In  parameterizing  the  JCZ3  equation  of  state,  one  must  find  a best 

attraction  parameter  (m)  and  a best  repulsion  parameter  (p)  for  the  equa- 

* 

tion  as  a whole  as  well  as  a set  of  potentials  (e,r  ) for  each  product 


species  encountered  in  the  calculation.  In  theory  it  should  be  possible 
to  determine  these  parameters  frcm  experimental  Hugoniot  data. 


r 

DISCUSSION 

L Because  of  the  limited  amount  of  data  used,  Mader's  parameterization1^ 

of  BKW  was  severely  restricted.  With  an  expanded  data  base  a more  severe 

test  of  the  parameterization  will  result.  Since  Mader's  parameterization 

of  BKW  in  1963,  much  useful  detonation  data  has  appeared  in  literature. 

Not  only  do  detonation  pressure  and  detonation  velocity  data  for  a large 

number  of  explosives  abound  in  the  literature1^' but  in  recent  years 

18  19  20  21 

several  detonation  product  composition  studies  ' ' ' have  also 

appeared.  The  approach  used  was  to  select  several  explosives  displaying 
a wide  variety  of  elemental  compositions  on  which  experimental  detonation 
pressure  as  a function  of  loading  density,  detonation  velocity  as  a 
function  of  loading  density,  and  detonation  product  composition  data  was 
available.  The  explosives  selected  comprised  both  C,H,N,0  and  non- 
C,H,N,0  explosives.  In  Table  III  are  listed  the  explosives  used  in  this 
study  along  with  their  elemental  compositions  and  experimental  data  avail- 
able in  the  literature.  This  set  of  explosives  should  comprise  a very 


severe  test  for  any  set  of  BKW  parameters  employed. 

When  the  experimental  detonation  product  composition  and  the  BKW 

calculated  Chapman-Jcuget  (GJ)  i sen tropes  are  compared  for  both  HMX  and 

18 

TNT,  an  interesting  observation  arises.  Omellas  has  noted  that  the 
detonation  products  resulting  from  heavily  confined  charges  of  HMX  and 
TOT  attain  equilibrium  under  nonideal  gas  conditions'  between  approxi- 
mately 1500°K  and  1800°K.  The  RUBY  code  was  used  to  calculate  the  CJ 
i sen  tropes  at  these  two  temperatures  for  HMX  and  TOT.  The  experimental 
product  composition  should  then  fall  within  the  range  of  values  for  the 


T&ble  Ilia.  Explosives  Used  in  the  Reparameterization  of  the  BKW  Equation  of  State  Along  With 


1500°K  and  1800°K  calculated  i sen  tropes  since  in  this  range  the  product 
species  composition  is  "frozen".  The  comparison  can  be  seen  in  Table  IV 
for  HMX  and  TNT.  In  both  cases  it  is  iimediately  obvious  that  the  calcu- 
lated quantity  of  CO  is  much  smaller  than  the  observed  quantity  while  the 
calculated  amount  of  C02  is  much  larger  than  the  observed  amount  of  this 
product  species.  Furthermore,  the  observed  amount  of  falls  very 
near  the  lower  end  of  the  calculated  range  for  this  product  species. 

The  relative  quantities  of  these  three  carbon-containing  species  are 
governed  by  equilibrium 

200  (g)" >002(g)  +C(s)  (5) 

Mader  assumed  formation  of  carbon  in  the  standard  state  (i.e. , AHf  = 0) ; 
however,  it  may  be  noted  that  the  above  equilibrium  can  be  shifted  to  the 
left  by  increasing  the  heat  of  formation  of  C^-  which  is  input  into  the 
computer.  The  BKW  calculated  i sen  tropes  would  then  contain  more  CO  and 
less  002  and  C(s)*  This  is  in  the  correct  direction  to  satisfy  the 
deficiencies  discussed  above.  It  was  thus  decided  to  add  an  additional 
parameter  - the  heat  of  formation  of  C^  - in  the  reparameterization  of 
the  BKW  equation  of  state.  Added  impetus  to  this  decision  was  the  obser- 
vation that  the  carbon  produced  was  not  graphite  but  an  amphorous  form. 
Additionally,  several  other  researchers  (particularly  scientists  at 
Lawrence  Livermore  laboratory  and  in  the  Soviet  union)  were  using  positive 
heats  of  formation  for  carbon  in  their  computer  codes. 

There  are  very  few  measurements  of  detonation  temperatures  at  the 
CJ  state  and  these  have  relatively  large  uncertainties^ ' ^ . Nevertheless, 
it  was  intuitively  felt  that  the  calculated  detonation  tenperatures  should 
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be  higher.  For  a fixed  set  of  the  other  BKW  parameters,  (a,  0,  <,  k^'s), 
the  larger  the  value  of  the  parameter  0,  the  higher  the  calculated 
detonation  tsnperature. 

The  approach  taken,  therefore,  in  reparameterizing  the  BKW  equation 
of  state  was  the  systematic  variation  of  the  five  enp.irical  parameters 
(a,  8,  k,  9 and  AH^)  while  providing  for  the  dependence  of  detonation 
temperature  on  0 and  the  dependence  of  equation  (5)  on  the  heat  of  -forma- 
tion of  C^sy.  Additionally,  each  product  species  has  a characteristic 
covolume,  k^.  These  product  species'  covolumes  were  also  treated  as 
empirical  parameters;  however,  they  were  constrained  within  certain, 
reasonable  limits  and  not  allowed  to  completely  vary.  In  the  parameteri- 
zation, not  only  were  detonation  pressures  and  velocities  used  in  the 
fit,  but  also  product  composition.  Furthermore,  the  fit  was  extended  to 
include  the  non-C,H,N,0  explosives  BTF,  KX-23M,  RX-23AB  and  RX-23AC  as 
well  as  the  C,H,N,0  explosives  HMX,  TNT  and  BTNEA.  These  extensions 
should  give  a very  severe  test  of  the  BKW  equation  of  state  both  from 
the  standpoint  of  the  wide  range  of  explosive  properties  used  in  the 
parameterization  as  well  as  the  types  of  explosives  employed  in  the  fit 
attaipt. 

The  results  of  this  "brute  force"  parameterization  attempt  are  given 
in  Tables  V and  VI.  The  "beet"  fit  set  of  parameters  and  product  specie 
colvolunes  are  displayed  in  Table  VII.  The  results  shown  in  Tables  V and 
VI  are  exceptionally  good  considering  the  wide  range  of  detonation  proper- 
ties fit  and  the  various  types  of  explosives  used.  Comparison  of  calcu- 
lated results  with  experimental  data  danonstrates  that  the  differences  be- 


tween the  two  are  usually  within  the  known  experimental  error  which  is 
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generally  no  better  than  5%.  Again  in  accordance  with  Omellas ' find- 
ings (18),  a "freeze-out"  temperature  range  of  1500°K-1800°K  has  been 
employed  for  calculated  detonation  product  composition  data. 

Because  the  JCZ3  equation  is  developed  from  first  principles  and, 
therefore,  the  parameters  do  have  physical  significance,  it  was  felt  that 
a purely  "brute  force"  type  of  fit  would  be  unsatisfactory.  The  values 
assigned  to  the  parameters  should  be  determined  from  experimental  data 
and,  therefore,  should  be  physically  meaningful.  Fortunately,  for  many 
of  the  detonation  product  species,  experimental  Hugoniot  data  does  exist. 
For  each  of  these  product  species  then,  the  pair  potential  parameters 
(e/k,  r*)  can  be  determined.  First,  however,  the  attraction  and  repul- 
sion parameters  (m,p)  of  the  JCZ  equation  of  state  must  be  set.  Because 
this  equation  of  state  was  ideally  written  for  spherical  molecules  in  a 
face-centered  cubic  (fee)  lattice  structure,  it  is  necessary  to  find 
experimental  Hugoniot  data  for  a spherical  molecule  with  a face-centered 
cubic  lattice  structure  that  can  be  used  to  set  the  value  of  the  attrac- 
tion and  repulsion  parameters  as  well  as  its  cwn  pair  potential  parameters. 
A molecule  for  which  much  experimental  data  is  available  and  which  sat- 
isfies the  above  requirements  is  argon  (24-28) . Once  the  attraction  and 
repulsion  parameters  are  determined  from  argon  Hugoniot  data,  they  are 
used  unchanged  in  the  determination  of  the  pair  potential  parameters  for 
each  of  the  detonation  product  species.  Because  the  product  species  are 
not  generally  of  a spherical  shape  and  often  do  not  pack  in  a face-centered 
cubic  lattice  structure,  stall  deviations  from  the  "true"  pair  potential 
values  are  expected.  The  detonation  product  species  for  which  experimental 
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Hugoniot  data  are  available  are  water  (H20)  (29),  nitrogen  (N2)  (30,31), 
carbon  dicocide  (OC>2)  (30),  oxygen  (C>2)  (32),  and  hydrogen  (H2)  (28,  33-35). 

The  procedure  used  in  establishing  the  values  of  the  attraction  and 
repulsion  parameters  and  the  pair  potential  parameters  was  to  use  the 
TIGER  computer  code  to  calculate  a Hugoniot  curve  for  the  species  of 
interest.  This  curve  was  compared  with  the  experimental  Hugoniot  curve. 
Adjustments  were  made  to  the  parameters  until  the  fit  was  within  experi- 
mental error.  By  working  with  Hugoniot  data  in  setting  the  parameter 
values  as  opposed  to  actual  detonation  property  data  of  explosives,  one 
restricts  the  number  of  variables  that  must  be  simultaneously  adjusted  to 
obtain  a reasonable  fit.  With  the  exception  of  argon,  only  two  variables 
must  be  adjusted  to  obtain  the  pair  potential  parameters  for  any  one  species. 
If  detonation  property  data  were  employed,  one  would  be  required  to  simul- 
taneously adjust  all  of  the  pair  potential  parameters  of  the  product 
species  as  well  as  the  attraction  and  repulsion  parameters.  This  vrould  be 
a very  laborious  task.  Furthermore,  it  is  very  likely  that  the  "best"  fit 
pair  potential  parameters  could  no  monger  be  related  to  physical  properties 
of  the  system  since  there  would  be  a tendency  to  spread  deviations  over  all 
parameters.  It  is  essential,  therefore,  that  care  be  taken  to  insure 
realistic  values  for  the  pair  potential  parameters.  Estimates  of  these 
parameters,  based  on  physical  measurements,  can  be  obtained  from  several 
sources  (13,36,37) . It  should  be  renarbered  that  the  calculated  parameters 
will  vary  slightly  from  the  measured  values  because  of  the  restrictions 
of  the  JCZ3  equation  of  state.  Once  the  pair  potential  parameters  of 
species  far  which  experimental  Hugoniot  data  exist  have  been  determined, 
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then  these  parameters  for  the  other  product  species  (CH4,  NH^,  NO,  00)  must 
be  evaluated.  This  can  be  accomplished  by  obtaining  initial  estimates 
(13,36,37)  of  the  parameters  and  then  using  detonation  product  data  from 
several  explosives.  As  new  product  species  are  added  to  the  TIGER  library, 
it  will  be  necessary  to  determine  the  pair  potential  parameters. 

In  order  to  facilitate  the  parameterization  of  the  JCZ3  equation  of 
state,  a subroutine  was  incorporated  into  the  TIGER  computer  program. 

This  subroutine,  based  on  the  method  of  Bevington  (38),  simultaneously 
varies  either  two  (e/k>  r*,  p) , three  (e/k»  r* , p) , or  four  parameters 
(e/k,  r*,  p,  m)  until  a "best"  fit  set  of  parameters  are  obtained.  The 
"best"  fit  parameters  are  ascertained  by  a comparison  of  the  experimental 
Hugoniot  data,  which  is  input,  and  the  calculated  Hugoniot  data.  A 
"quality  of  fit"  parameter  was  used  to  determine  when  a satisfactory  fit  of 
the  Hugoniot  data  is  obtained.  For  all  species  except  argon  only  the  two- 
parameter  option  was  arployed.  With  argon,  both  the  three-  and  the  four- 
parameter  options  were  used  in  an  attempt  to  obtain  a satisfactory  fit. 

When  the  three-parameter  option  was  employed,  the  value  of  the  attraction 
parameter  was  set  equal  to  6.0  for  which  there  appears  to  be  seme  theoreti- 
cal justification  (36) . 

This  attempt  to  parameterize  the  JCZ3  equation  of  state  has  demon- 
strated two  major  problan  areas.  In  the  analysis  of  argon  attempts  to 
vary  all  four  parameters  (c/k,  r*,  p,  m)  met  with  no  success  - primarily 
because  the  tanperature  of  the  system  reached  extremely  high  values  during 
the  calculations  causing  the  program  to  terminate  with  an  error.  This  is 
strictly  a software  problem  - not  a theoretical  one.  When  the  attraction 
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parameter  vras  set  at  6.0  and  the  other  three  parameters  for  argon  were 
varied,  a minimum  chi  square  value  was  obtained.  At  this  minimum  the 
repulsion  parameter  was  14.330,  eA  was  48.98  and  r*  was  3.93.  However, 
when  the  attraction  parameter  was  set  at  6.0,  the  value  of  the  repulsion 
parameter  was  readjusted  for  each  run  and  only  two  parameters  were  simul- 
taneously varied,  no  overall  minimum  could  be  located.  These  results  are 
shewn  in  T&ble  VIII.  The  value  of  the  chi  square  appears  to  vary  inversely 
with  the  value  of  the  repulsion  parameter.  At  values  of  the  repulsion 
parameter  above  14.5,  the  calculated  temperatures  become  extremly  high 
and  the  program  terminates  by  error.  Furthermore,  the  values  of  the  pair 
potential  parameters  become  very  unrealistic  for  high  values  of  "p" . No 
explanation  is  forthcoming  about  why  the  two  points  at  "p"  equal  to  14.33 
are  not  equivalent.  It  is  most  certainly  due  to  a peculiarity  in  the 
fitting  subroutine.  The  second  problem  area  concerns  the  fit  of  the  cal- 
culated Hugoniot  curve  to  the  experimental  Hugoniot  curve.  The  data  shewn 
in  Table  IX  is  typical  of  that  obtained  on  all  fit  attempts.  The  calcu- 
lated and  experimental  Hugoniot  curves  typically  fit  very  well  at  the  high 
pressure  segment  of  the  curve;  however, at  the  lower  pressures  (below  30 
kiloatmospheres)  the  fit  becomes  increasingly  unsatisfactory.  This  result 
demonstrates  the  need  to  reexamine  the  JCZ3  equation  of  state  particularly 
at  the  lew  pressure  end  where  most  of  the  detonation  products  cure  in  the 
gas  phase  or  between  the  gas  and  normal  liquid  state.  Because  no  values 
of  the  attraction  and  repulsion  parameters  could  be  obtained,  the  parame- 
terization attempt  is  suspended  until  the  reevaluation  of  the  JCZ3 
equation  of  state  is  accomplished. 


Table  VIII.  Attempted  Fit  for  Liquid  Argon  Hugoniot3  Using  Simultaneous 
Variation  of  e/k  and  r*  With  m = 6.0  and  Various  Input  Values 
of  p. 


p 

cA 

r* 

Chi  Square 

12.0 

450.0 

3.34 

26.51 

12.50 

315.1 

3.43 

21.28 

13.00 

210.3 

3.53 

17.91 

13.50 

136.5 

3.65 

15.73 

13.60 

125.0 

3.67 

15.41 

13.70 

114.1 

3.70 

15.10 

13.90 

94.2 

3.75 

14.57 

14.00 

86.4 

3.77 

14.33 

14.33 

63.2 

3.87 

13.70 

14.50 

53.5 

3.91 

13.43 

aReference  26. 


Table  IX.  Typical  Fit  Atterrpt  of  Calculated  and  Experimental  Hugoniot 
Curves  for  Liquid  Argon3  (eA  ■ 48.98,  r*  *=  3.93,  p = 14.33, 
m = 6.00) 


Input  Volune 
(cc/gm) 

Exp.  Pressure 
(atmos.) 

Cede.  Pressue 
(atmos) 

% Error 

0.5585 

6713. 

12253. 

45.2 

0.5259 

11727. 

16535. 

29.1 

0.5026,. 

17965. 

21141. 

15.0 

0.4851 

25426, 

25955. 

' 2.0 

0.4780 

30000. 

28385. 

-5.7 

0.4600 

40000. 

36189. 

-10.5 

0.4460 

50000. 

44545. 

-12.2 

0.4316 

59421. 

56210. 

-5.7 

0.4184 

74064. 

71090. 

-4.2 

0.4071 

90248. 

88163. 

-2.4 

0.3975 

107973. 

107400. 

-0.5 

0.3892 

127237. 

128681. 

1.1 

0.3819 

148042. 

151793. 

2.5 

0.3754 

170386. 

176426. 

3.4 

0.3697 

194272. 

202156. 

3.9 

0.3646 

219697. 

228480. 

3.8 

0.3600 

246663. 

254819. 

3.2 

0.3558 

275169. 

280618. 

1.9 

0.3520 

305215. 

305429. 

0.0 

0.3486 

336801. 

328938. 

-2.4 

0.3454 

369928. 

351005. 

-5.4 

aPefenence  26. 
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CONCLUSION 

The  BKW  equation  of  state  has  been  successfully  reparameterized . ihe 
reparameterized  version  of  this  equation  of  state  is  far  superior  to  any 
previous  versions  in  that  it  not  only  satisfactorily  predicts  detonation 
pressures  and  detonation  velocities  for  ideal  (C,H,N,0)  explosives,  but 
also  predicts  these  properties  as  well  as  detonation  product  compositions 
for  both  ideal  and  non-ideal  explosives.  The  two  major  changes  in  the 
reparameterization  were  a significant  increase  in  the  value  of  the 
parameter  0 and  the  addition  of  a fifth  parameter  - the  heat  of  formation 
of  amorphous  carbon.  In  considering  the  data  in  Table  V,  one  notes  that 
BTF  is  the  only  explosive  where  significant  deviations  from  the  experimental 
data  are  calculated.  It  is  believed  that  for  BTF  the  "freeze  out"  tenpera- 
ture  range  may  be  significantly  higher  than  the  1500°-1800°K  range  employed. 
BTF  appears  to  be  a much  hotter-burning  explosive  than  any  of  the  other 
explosives  used.  As  shewn  in  Table  V,  when  a 2150°-2550°K  "freeze-out" 
range  is  used,  the  fit  is  very  good.  Another  point  should  be  noted  regard- 
ing this  reparameterized  BKW  equation  of  state.  After  the  reparameteri- 
zation was  completed,  an  error  in  one  of  the  methane  thermodynamic 
constants  (15)  was  discovered.  The  thermodynamic  constants  for  all 
product  species  are  given  in  Appendix  A.  As  long  as  the  reparameterized 
BKW  equation  of  state  is  used  with  this  incorrect  value,  the  results 
should  be  satisfactory.  Hie  new  parameters  compensate  for  the  error  in 
the  incorrect  methane  thermodynamic  constant.  If  this  error  is  corrected, 
a new  parameterization  must  be  acccnplished . 
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Attarpts  to  parameterize  the  JCZ3  equation  of  state  were  unsuccessful. 
It  is  felt  that  the  primary  reason  for  this  is  that  the  equation  itself 
needs  to  be  revised  in  order  to  better  calculate  detonation  properties  at 
lcwer  pressures.  Once  the  equation  has  been  revised,  attempts  at  parame- 
terization of  the  equation  can  continue . Once  satisfactorily  parameterized, 
the  JCZ3  equation  of  state  should  be  far  superior  to  any  other  equation 
of  state  currently  being  used. 
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